ZEROS OF RANDOM POLYNOMIALS ON C m 

THOMAS BLOOM AND BERNARD SHIFFMAN 

Abstract. For a regular compact set K in C™ 1 and a measure fi on K satisfying the 
Bernstein-Markov inequality, we consider the ensemble Vn of polynomials of degree N, 
endowed with the Gaussian probability measure induced by L 2 (/x). We show that for large 
N, the simultaneous zeros of m polynomials in Vn tend to concentrate around the Silov 
boundary of K\ more precisely, their expected distribution is asymptotic to N m fi eq , where 
fx eq is the equilibrium measure of K. For the case where K is the unit ball, we give scaling 
asymptotics for the expected distribution of zeros as N — > oo. 



1. Introduction 

A classical result due to Hammersley |Haj (see also |SVj ). loosely stated, is that the zeros 
of a random complex polynomial 

TV 
j=0 

mostly tend towards the unit circle \z\ = 1 as the degree iV —>■ oo, when the coefficients Cj 
are independent complex Gaussian random variables of mean zero and variance one. In this 
paper, we will prove a multivariable result ( Theorem 13. lj) . a special case ( Example 13. 5j) of 
which shows, loosely stated, that the common zeros of m random complex polynomials in 

C m , 

/*(*)= E ck J z i-'- z t forfc = l,...,m, (2) 

\J\<N 

tend to concentrate on the product of the unit circles \zj\ = 1 (j — l,...,m) as N —>■ oo, 
when the coefficients c 1 } are i.i.d. complex Gaussian random variables. 

The following is our basic setting: We let if be a compact set in C m and let /i be a 
Borel probability measure on K. We assume that K is non-pluripolar and we let Vk be its 
pluricomplex Green function. We also assume that K is regular (i.e., Vk = V K ) and that 
H satisfies the Bernstein-Markov inequality (see $2}. We give the space Vn of holomorphic 
polynomials of degree < iV on C m the Gaussian probability measure induced by the 
Hermit ian inner product 

(/,<?)= / fgdfi. (3) 

J K 

The Gaussian measure 7/v can be described as follows: We write / = Yl'f^i c jPv wnere {Pj} 
is an orthonormal basis of Vn with respect to J5J) and d(N) = dimP/v = C^" 1 ). Identifying 
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/ G Vn with c = (ci, . . . , Q(jv)) £ we have 

1 _| 12 

rf 7iv(s) = ^ye |c| dc. 

(The measure 7jv is independent of the choice of orthonormal basis {pj}-) In other words, a 
random polynomial in the ensemble (Vn, In) is a polynomial / = ^ . Cj p^-, where the Cj are 
independent complex Gaussian random variables with mean and variance 1. 

Our main result, Theorem 13.11 gives asymptotics for the expected zero current of k i.i.d. 
random polynomials (1 < k < to). In particular, the expected distribution E(Z/ lr ..j TO ) of 
simultaneous zeros of m independent random polynomials in (Vn,7n) has the asymptotics 

^ E (Z/i,...,/J ^Veq weak*, (4) 

where = (^ddVx) m is the equilibrium measure of K. Here, E(X) denotes the expected 
value of a random variable X. 

The reader may notice from (J3J) that the distributions of zeros for the measures on Vn 
considered here are quite different from those of the SU(m + 1) ensembles studied, for exam- 
ple, in |SZ1) ISZ4) IBSZ1) IBSZ2| IDS) . The Gaussian measure on the SU(m + 1) polynomials 
is based on the inner product 



(/) 9)n — I F N Gn, 

J5 2m + 1 



where F N , Gn € C[z , z 1} . . . , z m ] denote the degree iV homogenizations of / and g respec- 
tively. It follows easily from the SU(m + l)-invariance of the inner product that the expected 
distribution of simultaneous zeros equals ^-uj 171 (exactly), where u is the Fubini-Study 
Kahler form (on C m C CP m ). We note that, unlike (JHJ), this inner product depends on 
N; indeed, \\z J \\% = mi{N ^ J ^" iml |SZH (30)]. 

In this paper, we also give scaling limits for the expected zero density in the case of the 
unit ball in C m (Theorem 14.1)1 . The problem of finding scaling limits for more general sets 
in C m remains open. Another open problem is to establish the multivariable version of the 
following one variable result: For a regular subset K C C, it is known (see [SZ1, B12J) that 
with probability one, a sequence {/jv}jv=i,2,... °f random polynomials of increasing degree 
satisfies: 

lim —Z fN = fi eq weak* . 
n^oo yv 

2. Background 

We let L denote the Lelong class of plurisubharmonic (PSH) functions on C m of at most 
logarithmic growth at oo. That is 

£:={«£ PSH(C m ) | u{z) < log + + 0(1)} (5) 

For K a compact subset of C m , we define its pluricomplex Green function Vk{z) via 

Vk(z) = sup{u(z) | u G C, u < on K}. (6) 

We will assume K is regular, that is by definition, Vk is continuous on C m (and so Vk = V£, 
its uppersemicontinuous regularization). The function Vk is a locally bounded PSH function 
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on C m and, in fact 

V^-log + ||z|| =0(1) . (7) 

By a basic result of Bedford and Taylor [BTlJ (see K1J), the complex Monge- Ampere operator 
(dd c ) m = (2idd) m is defined on any locally bounded PSH function C m and in particular on 
Vk- The equilibrium measure of K is defined by (see jKTJ Cor. 5.5.3]) 

(\ m 
% -ddV K \ (8) 

Since Vk satisfies (J7|), it is a positive Borel measure, here normalized to have mass 1. The 
support of the measure fj, eq (K) is the Silov boundary of K for the algebra of entire analytic 
functions p3T2jj . In one variable, i.e. K C C, Vk is the Green function of the unbounded 
component of C \ K with a logarithmic pole at oo, and fi eq (K) = ^AVk, where A is the 
Laplacian |Raj . 

Let /i be a finite positive Borel measure on K. The measure /j, is said to satisfy a Bernstein- 
Markov (BM) inequality, if, for each e > there is a constant C = C(e) > such that 

\\p\\K<Ce £de ^\\p\\ LHp) (9) 

for all holomorphic polynomials p. Essentially, the BM inequality says that the L 2 norms and 
the sup norms of a sequence of holomorphic polynomials of increasing degrees are "asymp- 
totically equivalent" . 

The question arises as to which measures actually satisfy the BM inequality. It is a result 
of Nguyen- Zeriahi NZ] combined with jKll Cor. 5.6.7] that for K regular, \x eq {K) satisfies 
BM. This fact is used in Examples 13. 5fJ3~o1 In [Bll, Theorem 2.2], a "mass-density" condition 
for a measure to satisfy BM was given. (See also |BLj .) 

Our proof uses the probabilistic Poincare-Lelong formula for the zeros of random functions 
( Proposition 12. II below) . Considering a slightly more general situation, we let g\, . . . ,gd be 
holomorphic functions with no common zeros on a domain U C C m . (We are interested in 
the case where U = C m and {gj} is an ortho normal basis of Vn with respect to the inner 
product as discussed above.) We let T denote the ensemble of random holomorphic 
functions of the form f = Yl c j 9ji where the Cj are independent complex Gaussian random 
variables with mean and variance 1. We consider the Szego kernel 

d 

S T {z,w) = ^gj{z)gj(w) . 
i=i 

For the case where the gj are orthonormal functions with respect to an inner product on 
0(U), Sj?(z,w) is the kernel for the orthogonal projection onto the span of the gj. 

Under the assumption that the gj have no common zeros, it is easily shown using Sard's 
theorem (or a variation of Bertini's theorem) that for almost all fi, . . . , G J 7 , the differ- 
entials dfi, . . . , dfk are linearly independent at all points of the zero set 

loc(/i, ...,fk):= {z e U: h (z) = ■■■ = f k (z) = 0} . 

This condition implies that the complex hypersurfaces loc(/j) are smooth and intersect 
transversely, and hence loc(/i, . . . , /&) is a codimension k complex submanifold of XI . We 
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then let Zf lt ___j k G V' k,k (U) denote the current of integration over loc(/i, . . . , /&): 

f,-r) [ V, V eV m - k > m -\U). 

J\oc(f u ...,f k ) 

We shall use the following Poincare-Lelong formula from |SZ3|. ISZ4j: 

Proposition 2.1. The expected zero current ofk independent random functions fi, ■ ■ ■ , fk £ 
T is given by 

E(Z fu ...jJ=Q-dd\ogSAz,z?) ■ 

The proof follows by a verbatim repetition of the proof of Proposition 5.1 in SZ3J (which 
gives the case where the gj are normalized monomials with exponents in a Newton polytope). 
The codimension k — 1 case was given in {SZlJ (for sections of holomorphic line bundles), and 
in dimension 1 by Edelman-Kostlan |EKj . (The formula also holds for infinite-dimensional 
ensembles; see jSol ISZ4j .) The general case follows from the codimension 1 case together 
with the fact that 

E(Z fl _ fk ) = E(Z fl ) A • • • A E(Z fk ) = E(Z f ) k , (10) 

which is a consequence of the independence of the fj. The wedge product of currents is not 
always defined, but Z^ A ■ ■ • A Zf k is almost always defined (and equals Zf ly _j k whenever the 
hypersurfaces loc(/j) are smooth and intersect transversely), and a short argument given in 
[SZ3J yields (fTUjl . (In fact, the left equality of ([111)) holds for independent non-identically- 
distributed fj, as proven in |SZ3j .) We note that the expectations in (jlOj) are smooth forms. 

3. Random polynomials on polynomially convex sets 

Theorem 3.1. Let /i be a Borel probability measure on a regular compact set K C C m , 
and suppose that (K, fi) satisfies the Bernstein- Markov inequality. Let 1 < k < m, and let 
iV%il%) denote the ensemble of k -tuples of i.i.d. Gaussian random polynomials of degree 
< N with the Gaussian measure d-j^ induced by L 2 (fi). Then 

]^fc E 7^(%.-./fe) -»• \Z 9 ^ V k\ weak*, as N -> oo , 

where Vk is the pluricomplex Green function of K with pole at infinity. 
To prove Theorem 13.11 we consider the Szego kernels 

d(N) 

S N (z,w) := S {Vni1n) (z,w) = ^2p j (z)p j (w) , 

3=1 

where {pj} is an L 2 ( / u)-orthonormal basis for Vn- Our proof is based on approximating the 
extremal function Vk by the (normalized) logarithms of the Szego kernels Sn(z, z) (Lemma 

EM- 

We begin by considering the polynomial suprema 

*%{z) = bu P {|/(z)| :fEV N , \\f\\ K < 1} • (11) 
Since ^log/ G C, for / G Vn, it is clear that ^log$^ < Vk, for all N. Pioneering work 
of Zaharjuta |Zaj and Siciak |Sil| ISi2j established the convergence of -^log$^ to Vk- The 
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uniform convergence when K is regular seems not to have been explicitly stated and we give 
the proof below. 

Lemma 3.2. Let K be a regular compact set in C m . Then 

iiog$£(s)-yjc(z) 

uniformly on compact subsets ofC m . 

Proof. We first note that 1 < $j < < &j+k, for j, k > 0. By a result of Siciak |Silj 

and Zaharjuta |Zaj (see (KlJ Theorem 5.1.7]), 

V K (z) = lim ^log^(z) = supllog**(*) , (12) 

for all 2 G C m . 

We use the regularity of K to show that the convergence is uniform: let 

^ = llog$£>0. 

Thus for N, k > 1, j > 0, we have 

IpNk + j tpj < (Nk + j) ifjNk+j ■ 

Since ip^ < ipNk, we then obtain the inequality 

Nk j Nk , 

^ " NkT]^ + NkTj^ - ~Nk+~j ■ (13) 
Fix e > 0. For each a G C m , we choose iV a G Z + such that 

Vk-(a) - ipN a (a) < £ and ~TT~ < 5 ' 
and then choose a neighborhood U a of a such that 

\Vk(z) - V K (a)\ < e, ip Na (z) > i>N a (a) - e, V *lf < e, for z e U a . 

Now let iV > N*, and write iV = N a k + j, where k > N a , < j < N a . By f|T2]l - (Pjl . we 
have 

N k N 1 

0<V K -^ N <V K - W f—.^N a <V K - -rr^^ <V K - ^ Na + -tt—Vk . (14) 

N a k + j N a + 1 N a + 1 

Hence, for all iV > N% and for all z G U a , we have 

< V K {z) - i> N (z) < V K {z) - 1> Na (z) + e 

= [V K {a) - ip Na (a)] + [V K {z) - V K {a)\ + \^ Na (a) - iJj Na (z)} + e 
< As. (15) 
Hence for each compact A C C m , we can cover A with finitely many U ai , so that we have by 

\\V K - if> N \\ A < 4e V iV > maxJV?. . 

i 1 

□ 
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Lemma 3.3. For all e > 0, there exists C = C £ > such that 

1 <|^|<Ce^(iV). 



d(N) ~ $Uz) 2 
Proof. Let f e V N with \\f\\ K < 1. Then 



Sn(z, w)f(w) dfj,(w) 



K 



< / \S N (z,w)\d/j,(w) 



K 

i 



< / S N (z,z)2S N (w,w)2 dfi(w) = S N (z,z)^ \\S N (w,w)2\\ L i {tl) 



K 



< S N (z,z)z \\l\\ L 2 M \\S N (w,w)2\\ L 2 M = S N (z,z)i d(N)» . 

Taking the supremum over / G Vn with \\f\\K < 1, we obtain the left inequality of the 
lemma. 

To verify the right inequality, we let {pj} be a sequence of L 2 (/i)-orthonormal polynomials, 
obtained by applying Gram-Schmid to a sequence of monomials of non- decreasing degree, 
so that {pi, . . -PdiN)} is an orthonormal basis of Vn (for each N e Z + ). By the Bernstein- 
Markov inequality Q, we have 

\\ Pj \\ K < Ce £deg ^ 

and hence 

\p 3 {z)\ < WptU^^z) < Ce £de ^ &< cgpj (z) < Ce* N $%{z) , for j < d(N). 
Therefore, 

d(N) 

S N (z,z) = \PM 2 < d{N)C 2 e 2F ' N ^ K N {zf . 

Lemma 3.4. Under the hypotheses of Theorem YS. 1\ we have 

log S N {z,z) -> V K {z) 

uniformly on compact subsets ofC m . 

Proof. Let e > be arbitrary. Recalling that d(N) = ( N ^ n ) } we have by Lemma 13.31 

m , ,, T x 1, f Sn(z,z)\ logC m, ,, T . 

-- log(iV + m) < - log < + s + - log(JV + m) . 

Since e > is arbitrary, we then have 

1 / S N (z,z) 



□ 



(16) 



The conclusion follows from Lemma 13.21 and (fTB|) . □ 

Proof of Theorem 1.9. it It follows from Lemma 13.41 and the fact that the complex Monge- 
Ampere operator is continuous under uniform limits [BT1 , 

\ k / • \ k 



27rA^ ° ' ' ' 7 V 7T 



ddlog Sn(z, z) — ► — 99Vx(^) I u>ea/c 
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3 I ' 



The conclusion then follows from Proposition 12.11 

Example 3.5. Let K be the unit polydisk in C m . Then V K 
boundary of K is the product of the circles \zj\ = 1 (J = 1, 
where dOj is the angular measure on the circle \zj\ = 1. 

The monomials z J := z{ x ■■■z J r £ l , for |J| < iV, form an orthonormal basis for V 
random polynomial in the ensemble is of the form 



□ 

the Silov 



, m 



max™^ log 

and dfi eq = {j^) m dQi ■ ■ ■ d9. 



N ■ 



A 



\J\<N 

where the cj are independent complex Gaussian random variables of mean zero and variance 
one. By Theorem 13.11 E 7 & (Zf lt _j m ) — > (^) m d9i ■ ■ ■ d6 m weak*, as — ► oo. In particular, 
the common zeros of m random polynomials tend to the product of the unit circles \zj\ — 1 
for j = 1, 



m. 



Example 3.6. Let K be the unit ball {\\z\\ < 1} in C m . Then the Silov boundary of K is 
its topological boundary {\\z\\ = 1}, Vk(z) = log" 1 " \\z\\, and \i eq is the invariant hypersurface 
measure on II z II = 1 normalized to have total mass one. 



4. Scaling limit zero density for orthogonal polynomials on S 2m ~ x 

Examples 13.51 and 13.61 both reduce to the unit disk in the one variable case. In that case, 
detailed scaling limits are known (see, for example, IZJ) . For a more general compact set 
K C C with an analytic boundary, scaling limits are found in |SZ2j . 

In this section, we consider the case where K = {z 6 C m : \\z\\ < 1} is the unit ball and 
fi is its equilibrium measure, i.e. invariant measure on the unit sphere S 2 ™' 1 . We have the 
following scaling asymptotics for the expected distribution of zeros of m random polynomials 
orthonormalized on the sphere: 

Theorem 4.1. Let {V^,j^}) denote the ensemble of m-tuples of i.i.d. Gaussian random 
polynomials of degree < N with the Gaussian measure djjy induced by L 2 (S' 2m_1 , [f), where 
fi is the invariant measure on the unit sphere S 2m ~ l C C m . Then 



D N (log 



'Ml 



where 
1 



u 



, m— 1 



log 



d 



m—1 



du 



m—1 



e u - 1 



Proof. We write 



z — z 



ji 



An easy computation yields 

/ \z J \ 2 dfi(z) 
Js 2 ™- 1 

where 



Z = (z\, . . .,Zm), J = (jl, • • .,j m ) 

(m-l)\j 1 \--- Jm \ 1 



\J\ = h + '-- + 3 



(\J\+m-l)\ ('tTW) ' 
|J|! 



(17) 
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Thus an orthonormal basis for Vn on S 2 " 1 ^ 1 is: 



i 1. 

2 



We have 

N 



|J|<V k=0 v 7 |J|=fc v 7 

" //„ , ™ 1 \ 

|2fc 



y~v / A' - - /// — 1 



A:=0 



m — 1 



Hence 



SV(z, z) = g N (\\z\\ 2 ) , where g N (x) = ^T ( k + m ] - k 



x k . (19) 



k=0 

We note that 

i(m— 1) 



1 ^ m _n , „ , , l-x N+m 



9n = 7 7T7 G ( N , where G N (x) = — 

(m — 1)! 1 — 

We denote by O(j^) any function X(N, u) = \n{u) : Z + x IR ^ M satisfying: 



Vi? > 0, Vj G N, 3Cjy G M+ such that sup \\%\u)\ < ^ . (20) 



\u\<R 



We note that 



and hence 



Thus we have 



Wlog(l + ^)=» + »'o(l) 



(for |-u| < A/") , 



( 1 + l)" = e " + " 2o/1 



AT 



Hence 



AT 

Therefore 

7m- 1)! 



W 1+ *)-^ +0 (*)- (21) 

1 u\ 1 d m ~ x fe u -l\ „fl\ 

^ 9 H 1+ iv) = (^^^(— J +0 UJ • (22) 



log 



F m (u) + O ( ± ) , (23) 



AT" 

where F m is given in the statement of the theorem. 

Since the zero distribution is invariant under the SO(2m)-action on C m , we can write 



\dd\\z\\^ . (24) 
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Then D N (j^) is the density at the point 



. N 



u/2N 



m 



lz Nu2 = e u/N 



We shall compute using the local coordinates Q = pj + i6j = log Zj. Let 



2 



By Proposition 12. II and (|19|). we have 



—Yde ( l ^ 
2tt J \2 dpjdp k 



We note that 



n = m r 



1 ^ 



at the point z 



N 



(25) 



(26) 



We let 1 denote the m x m matrix all of whose entries are equal to 1 (and we let / denote 
the m x m identity matrix). By (|23|) and (|25J ) - (J2()j) . we have 

1 N " 



D 



N 



2^ 



1 + 



N 



Therefore, 



7T 



x det (2m~ 2 e 2u/N (\ogg N )"(e u/N ) 1 + 2 m" 1 e u/N (log g N )'(e u/N ) I 
det (m- l N 2 F'^u) 1 + N F' m (u) A . 



1 + ° { N 



D N (-) = 



N 



1 



1 + 0{ N 



+ m[m- l N 2 F'^u)] [NF^ 1 ^ 



= ±F^u)F^ur-i + o(± 



□ 



Remark: There is a similarity between the scaling asymptotics of Theorem 14.11 and that 
of the one-dimensional SU(1,1) ensembles in |BRj with the norms = ( L_1+ ' ? ) ^ 2 , for 
L G Z + . Then the expected distribution of zeros of random SU(1, 1) polynomials of degree 
N has the asymptotics fBRj Th. 2.1]: 



E N (Z f ) = D N (log \z\ 2 ) -dz A dz , 



where (in our notation) 



u 



10 



THOMAS BLOOM AND BERNARD SHIFFMAN 



References 

[BTl] E. Bedford and B. A. Taylor, A new capacity for plurisubharmonic functions, Acta Math. 149 
(1982), 1-40. 

[BT2] E. Bedford and B. A. Taylor, Fine topology, Silov boundary and (dd c ) n , J. Fund. Anal. 72 (1987), 
225-251. 

[BR] P. Bleher and R. Ridzal, SU(1, 1) random polynomials, J. Statist. Phys. 106 (2002), 147-171. 
[BSZ1] P. Bleher, B. Shiffman and S. Zelditch, Poincare-Lelong approach to universality and scaling of 

correlations between zeros, Comm. Math. Phys. 208 (2000), 771-785. 
[BSZ2] P. Bleher, B. Shiffman and S. Zelditch, Universality and scaling of correlations between zeros on 

complex manifolds, Invent. Math. 142 (2000), 351-395. 
[Bll] T. Bloom, Orthogonal polynomials in C", Indiana Univ. Math. J. 46 (1997), 427-452. 
[B12] T. Bloom, Random polynomials and Green functions, Int. Math. Res. Not. 2005 (2005), 1689-1708. 
[BL] T. Bloom and N. Levenberg, Capacity convergence results and applications to a Bernstein-Markov 

inequality, Trans. Amer. Math. Soc. 351 (1999), 4753-4767. 
[DS] T.-C. Dinh and N. Sibony, Distribution des valeurs de transformations meromorphes et applications, 

Comment. Math. Helv. 81 (2006), 221-258. 
[EK] A. Edelman and E. Kostlan, How many zeros of a random polynomial are real? Bull. Amer. Math. 

Soc. 32 (1995), 1-37. 

[Ha] J. H. Hammcrsley, The zeros of a random polynomial. Proceedings of the Third Berkeley sympo- 
sium on Mathematical Statistics and Probability, 1954-55, vol II, University of California Press, 
California 1956, pp 89-111. 

[IZ] I. Ibragimov and O. Zeitouni, On roots of random polynomials, Trans. Amer. Math. Soc, 349 
(1997), 2427-2441. 

[Kl] M. Klimek, Pluripotential Theory, London Math. Soc. Monographs, New Series 6, Oxford University 
Press, New York, 1991. 

[NZ] T. V. Nguyen and A. Zeriahi, Famine de polynomes presque partout bornees, Bull. Sci. Math.. 
Paris, 107 (1983), 81-91. 

[Ra] T. Ransford, Potential Theory in the Complex Plane, London Mathematical Society Student Texts 

28, Cambridge University Press, 1995. 
[SV] L. A. Shepp and R. J. Vanderbei, The Complex zeros of random polynomials, Trans. Amer. Math. 

Soc. 347 (1995), 4365-4384. 
[SZ1] B. Shiffman and S. Zelditch, Distribution of zeros of random and quantum chaotic sections of 

positive line bundles, Commun. Math. Phys. 200 (1999), 661-683. 
[SZ2] B. Shiffman and S. Zelditch, Equilibrium distribution of zeros of random polynomials, Int. Math. 

Res. Not. 2003 (2003), 25-49. 
[SZ3] B. Shiffman and S. Zelditch, Random polynomials with prescribed Newton polytope, J. Amer. 

Math. Soc. 17 (2004), 49-108. 

[SZ4] B. Shiffman and S. Zelditch, Number variance of random zeros, preprint 2005,|math.CV/0512652 



[Sil] J. Siciak, Extremal plurisubharmonic functions in C", Ann. Polon. Math. 39 (1981), 175-211. 

[Si2] J. Siciak, Extremal Plurisubharmonic Functions and Capacities in C™, Sophia Kokyuroku in Math- 
ematics, No. 14, Sophia Univ., Tokyo, 1982. 

[So] M. Sodin, Zeros of Gaussian analytic functions, Math. Res. Lett. 7 (2000), 371-381. 

[Za] V. P. Zaharjuta, Extremal plurisubharmonic functions, orthogonal polynomials, and the Bernstein- 
Walsh theorem for functions of several complex variables (Russian), Ann. Polon. Math. 33 
(1976/77), 137-148. 

Department of Mathematics, University of Toronto, Toronto, ON, CANADA M5S 3G3 
E-mail address: bloom@math.toronto.edu 



Department of Mathematics, Johns Hopkins University, Baltimore, MD 21218, USA 
E-mail address: shiffman@math.jhu.edu 



